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Gas  and  water  management  is  the  key  to  achieving  good  performance  from  a  polymer  electrolyte  mem¬ 
brane  fuel  cell  (PEMFC)  stack.  Imbalance  between  production  and  evaporation  rates  can  result  in  either 
flooding  of  the  electrodes  or  membrane  dehydration,  both  of  which  severely  limit  fuel  cell  performance. 
In  the  present  study,  a  mathematical  model  was  developed  to  evaluate  moisture  profiles  of  hydrogen 
and  air  flows  in  the  flow  field  channels  of  both  the  anode  and  the  cathode.  For  model  validation,  a  sin¬ 
gle  fuel  cell  was  designed  with  an  active  area  of  200  cm2.  Six  humidity  sensors  were  installed  in  the 
flow  fields  of  both  the  anode  and  the  cathode  at  457  mm,  1266  mm  and  2532  mm  from  the  inlets.  The 
experiment  was  performed  using  an  Arbin  Fuel  Cell  Test  Station.  The  temperature  was  varied  (25  °C, 
40  °C,  50  °C  and  60  °C),  while  hydrogen  and  air  velocities  were  fixed  at  3Lmin_1  and  6Lmin_1,  respec¬ 
tively,  during  the  operation  of  the  single  cell.  The  feed  relative  humidity  at  the  anode  was  fixed  at 
1.0,  while  the  feed  relative  humidity  at  the  cathode  was  fixed  at  0.005  (dry  air).  All  humidity  sensor 
readings  were  taken  at  steady  state  after  2  h  of  operation.  Model  predictions  were  then  compared  with 
experimental  results  by  using  the  least  squares  algorithm.  The  moisture  content  was  found  to  decrease 
along  the  flow  field  at  the  anode,  but  to  increase  at  the  cathode.  The  moisture  content  profile  at  the 
anode  was  shown  to  depend  on  the  moisture  Peclet  number,  which  decreased  with  temperature.  On 
the  other  hand,  the  moisture  profile  at  the  cathode  was  shown  to  depend  on  both  the  Peclet  number 
and  the  Damkohler  number.  The  trend  of  the  Peclet  number  in  the  cathode  followed  closely  that  of 
the  anode.  The  Damkohler  number  decreased  with  temperature,  indicating  increasing  moisture  mass 
transfer  with  temperature.  The  moisture  profile  models  were  successfully  validated  by  the  published 
data  of  the  estimated  overall  mass  transfer  coefficient  and  moisture  effective  diffusivity  of  the  same 
order  of  magnitude.  The  strategy  of  saturating  the  hydrogen  feed  and  using  dry  air,  as  in  the  present 
work,  has  been  shown  to  successfully  prevent  water  droplet  formation  in  the  cathode,  and  hence  prevent 
flooding. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Water  plays  a  critical  role  in  PEM  fuel  cells  (PEMFC).  Good  PEMFC 
performance  depends  on  good  water  management.  Extra  humidi¬ 
fication  of  the  reactant  gases  is  essential  in  PEMFC  operating  above 
60  °C  because  water  vapour  is  continuously  lost  in  the  exhaust 
gases,  and  membranes  must  be  hydrated  to  better  conduct  protons 
[1].  If  the  membrane  is  not  properly  hydrated,  its  ionic  resistance 
increases,  and  in  extreme  cases,  it  can  be  completely  destroyed. 
The  methods  normally  used  to  keep  the  membrane  hydrated  are 
internal  and  external  humidification.  In  “self-humidifying”  inter¬ 
nal  humidification,  a  catalyst  in  the  membrane  electrode  assembly 
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(MEA)  is  modified  to  not  only  to  retain,  but  also  to  produce  water 
[2].  This  process  is  passive  and  most  suitable  for  small  cells.  In  con¬ 
trast,  in  the  external  humidification  method,  the  fuel  gas  and  air  are 
passed  through  water  columns  of  separate  humidifier  bottles  until 
they  are  fully  saturated  before  entering  the  fuel  cell  stack  [3].  The 
humidifier  bottle  temperatures  are  controlled  independently  from 
the  stack  temperature  to  get  the  desired  gas  temperature  and  rela¬ 
tive  humidity.  The  external  humidification  method  is  widely  used  in 
small  scale  laboratory  fuel  cell  experiments  as  well  as  in  larger  pilot 
plant  scale  fuel  cell  systems  due  to  its  simplicity.  This  laboratory  has 
investigated  the  humidification  load  of  a  similar  external  humidi¬ 
fier  [4].  In  the  present  work,  a  fuel  cell  operating  technique  was 
developed  to  manage  water  in  the  fuel  cell  by  maintaining  a  mois¬ 
ture  gradient  across  the  membrane  from  the  anode  to  the  cathode. 
The  effectiveness  of  the  technique  was  investigated  via  modelling 
and  measuring  the  variation  of  humidity  in  the  air  and  hydrogen 
along  the  flow  fields  of  both  electrodes  in  a  single  cell  fuel  cell. 
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Nomenclature 

as  specific  surface  area  of  the  electrode  (m2  m-3 ) 

b  flow  channel  depth  (m) 

reference  molar  concentration  of  oxygen  (mol  m-3 ) 
dp  pore  diameter  (m) 

Da  Damkohler  number  (-) 

De  ff  diffusivity  (m  s-1 ) 

Do2  diffusion  coefficient  of  the  dissolved  oxygen  in  the 
electrolyte  (mol  cm-1  s-1) 

F  Faraday’s  constant  (A  s  mol-1 ) 

G  coefficient 

H  Henry’s  constant  (Pa  m3  mol-1 ) 

i’o  exchange  current  density  ( mA  cm-3 ) 

ioc  exchange  current  density  in  the  cathodic  active  layer 
(mA  cm-2) 

7a  moisture  condensation  flux  on  the  surface  of  the 

anode  (mol  cm-2  s-1) 

Jc  moisture  evaporation  flux  on  the  surface  of  the  cath¬ 

ode  (mol  cm-2  s-1) 

Jw  moisture  flux  across  the  membrane  (mol  cm-2  s-1 ) 

kA  moisture  mass  transfer  coefficient  at  the  anode 

(mol  cm-2  s-1) 

I<  K-v alue  of  moisture 

I<4  constant 

/<5  constant 

I<q  constant 

kc  moisture  mass  transfer  coefficient  at  the  cathode 

(mol  cm-2  s-1) 

kg  overall  moisture  mass  transfer  coefficient 

(mol  cm-2  s-1) 

L  total  length  of  the  flow  field  (m) 

m  evaporation  index  (-) 

M  coefficient  (-) 

n  evaporation  index  (-) 

p  total  pressure  (Pa) 

psat  saturated  vapour  pressure  of  moisture  (Pa) 

Q.  coefficient 

Pe  Peclet  number  (-) 

r  rate  of  water  production  by  the  electrochemical 

reaction  (mols-1) 

R  gas  constant  (J  mol-1  K-1 ) 

te  active  layer  thickness  (m) 

T  temperature  (K) 

v  gas  velocity  (ms-1) 

w  channel  width  (m) 

x*  equilibrium  moisture  content  of  membrane 

xci  molar  moisture  content  at  the  surface  of  cathode 

yA  molar  moisture  content  at  the  anode 

y*  equilibrium  moisture  content  at  membrane  surface 

yAF  hydrogen  feed  moisture  content 

yc  molar  moisture  content  at  the  cathode 

yCF  air  feed  moisture  content 

yci  molar  moisture  content  of  gas  at  the  surface  of  cath¬ 

ode 

z  distance  along  flow  field  (m) 

z*  dimensionless  distance 

Greek  letters 

£  dry  porosity  of  the  electrode 

0  relative  humidity 

A<peq  potential  difference  between  the  electrode  and  the 
membrane  (V) 

<Ps,  (pm  electrode  and  membrane  potential  (V) 


2.  Water  transport  in  PEMFC 

Membrane  hydration  is  affected  by  water  transport  in  the  mem¬ 
brane  itself,  which,  in  turn,  is  affected  by  the  condition  of  the  inlet 
gases  and  the  operating  parameters  of  the  fuel  cell.  Water  is  trans¬ 
ported  through  the  membrane  in  three  ways:  electro-osmotic  drag 
by  protons  from  the  anode  to  the  cathode,  back  diffusion  due  to  con¬ 
centration  gradients  from  the  cathode  to  the  anode  (or  vice  versa 
in  limited  cases),  and  convective  transfer  due  to  pressure  gradients 
within  the  stack. 

Current  research  has  focused  on  both  the  electro-osmotic  drag 
and  the  back  diffusion  as  the  dominant  mechanisms  for  water  trans¬ 
port  within  the  membrane.  Both  water  transport  mechanisms  are 
functions  of  the  fuel  cell  temperature,  current  density  and  mem¬ 
brane  water  content.  In  electro-osmotic  drag,  protons  that  move 
from  the  anode  to  the  cathode  drag  water  molecules  with  them. 
The  number  of  water  molecules  carried  by  each  proton  is  normally 
between  0.5  and  1.5  molecules  per  proton  [5].  As  a  result,  the  mois¬ 
ture  content  in  the  hydrogen  decreases  along  the  flow  channel. 
As  more  current  is  drawn  from  the  fuel  cell,  the  flow  of  protons 
and  water  from  the  anode  to  the  cathode  increases.  However,  at 
the  same  time,  water  that  is  being  produced  at  the  cathode  by  the 
oxygen  reduction  reaction  begins  to  diffuse  back  through  the  mem¬ 
brane  to  the  anode.  The  key  to  understanding  membrane  hydration 
is  the  balance  that  has  to  be  struck  between  the  electro-osmotic 
drag  and  back  diffusion  [5]. 

At  high  current  density,  the  membrane  dries  out  because  water 
transport  from  the  anode  by  electro-osmotic  drag  exceeds  its  trans¬ 
port  to  the  anode  by  back  diffusion  from  the  cathode  [6].  The  pores 
of  the  dry  membrane  also  shrink,  thus  further  limiting  the  back 
diffusion  of  water.  Water  transport  due  to  back  diffusion  is  not  suf¬ 
ficient  to  prevent  membrane  dehydration  [6,7].  The  net  water  flux 
from  the  anode  to  the  cathode  is  directly  proportional  to  the  current 
density  [7,8].  At  high  current  densities,  the  large  electro-osmotic 
drag  from  the  anode  dries  out  the  anode  and  floods  the  cathode 
[7,8].  The  reduction  reaction  at  the  cathode  is  adversely  affected  by 
the  increase  in  water  content  at  the  cathode  [8]. 

Cell  voltage  can  also  be  reduced  when  membrane  moisture  con¬ 
tent  is  low  due  to  anode  dehydration  [5].  High  cell  voltage  can  be 
achieved  by  saturating  the  inlet  hydrogen  with  water  to  increase 
the  water  flux  from  the  cathode  to  the  anode  [5].  Dry  hydrogen  car¬ 
ries  away  three  and  a  half  times  more  water  molecules  from  the 
anode  than  does  the  air  from  the  cathode  [9].  If  both  inlet  gases  are 
humidified,  the  water  profile  in  the  membrane  electrolyte  is  uni¬ 
form,  the  resistance  of  the  membrane  remains  low,  and  higher  cell 
voltages  can  be  achieved.  However,  the  saturated  air  at  the  cathode 
may  prevent  evaporation  of  the  water  produced  by  the  oxidation 
reaction  and  cause  flooding,  which  reduces  the  cell  voltage. 

Recent  works  on  water  management  have  shifted  their  focus 
of  managing  water  in  PEMFC  from  water  transport  in  the  MEA  to 
water  transport  in  the  gases,  which  can  be  controlled  by  manipulat¬ 
ing  operating  variables  of  the  PEMFC  stack  such  as  the  temperature, 
humidity  and  pressure  of  the  gases.  The  understanding  of  the  roles 
of  these  variables  on  water  management  is  still  rudimentary.  This 
is  reflected  in  the  use  of  either  dry  gases  or  fully  humidified  satu¬ 
rated  gases  [10-12].  The  water  content  of  the  gases  in  the  PEMFC 
cell  was  previously  found  to  affect  current  density  distribution  on 
the  electrode  and  therefore  the  PEMFC  performance  [13].  It  has  also 
been  suggested  that  effective  water  management  could  be  achieved 
by  controlling  water  content  in  the  feed  gases  [10].  However,  the 
low  current  density  and  performance  of  the  PEMFC  were  attributed 
to  low  proton  transport  caused  by  the  dryer  polymer  electrolyte 
membrane,  which  in  turn  was  caused  by  the  dryer  air  [10-12].  The 
conventional  solution  for  maintaining  a  sufficiently  wet  membrane 
for  maximum  proton  transport  is  to  use  fully  humidified  saturated 
hydrogen  and  air.  However,  this  solution  complicates  the  balance  of 
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plant  design  and  operation  and  often  leads  to  flooding  because  liq¬ 
uid  water  will  always  condense  in  the  air  at  the  saturated  condition 
approaching  water  dew  point. 

The  main  challenge  in  doing  this  is  maintaining  a  sufficiently 
wet  membrane  for  maximum  proton  transport  without  flooding  the 
PEMFC.  Flooding  occurs  when  the  wet  gas  in  the  PEMFC  reaches  its 
dew  point  [13]  or,  in  other  words,  its  saturation  point.  Humidifying 
both  gases  fully  to  saturation  will  certainly  lead  to  flooding.  In  the 
present  work,  an  alternative  operating  technique  for  better  water 
management  by  manipulating  the  water  content  and  temperature 
of  the  feed  gases  is  proposed.  Water  produced  at  the  cathode  is 
removed  as  soon  as  it  is  formed  through  evaporation  by  the  dryer 
air  to  avoid  flooding  at  the  cathode,  which  may  cause  the  MEA  to 
become  dry.  If  a  fully  humidified  and  saturated  hydrogen  is  used  at 
the  anode  and  a  partially  humidified  air  is  used  at  the  cathode,  the 
water  flux  is  forced  to  flow  from  the  anode  to  the  cathode  because 
of  the  difference  of  humidity  in  the  hydrogen  and  air  feed.  In  this 
manner,  the  water  content  in  the  membrane  can  be  maintained  and 
flooding  avoided. 


3.  Modelling  of  water  profile  in  the  PEMFC 

A  clearer  understanding  of  water  transport  in  a  PEM  fuel  cell  can 
be  obtained  by  validating  water  transport  models  with  experimen¬ 
tal  data.  Many  models  developed  to  explain  water  transport  in  PEM 
fuel  cells  in  the  past  were  either  based  on  water  flow  dominated  by 
transport  in  the  MEA  or  in  the  gas  flow  channels. 

Early  one  dimensional  models  considered  water  flow  across  the 
membrane  in  the  MEA  either  as  a  function  of  the  hydraulic  pres¬ 
sure  gradient  and  the  electro-osmotic  drag  in  a  uniformly  wetted 
membrane  with  porous  two  phase  catalyst  layers  [14]  or  as  a  func¬ 
tion  of  the  electro-osmotic  drag  and  the  diffusion  of  water  down  a 
moisture  content  gradient  across  the  electrolyte  membrane  [5].  The 
MEA  centred  models  were  varied  by  replacing  the  homogenous  two 
phase  model  of  catalyst  layer  by  an  agglomerate  model  [15]  and  by 
combining  the  electrolyte  membrane  model  and  the  agglomerate 
catalyst  layer  model  in  a  single  MEA  [16].  When  liquid  water  flood¬ 
ing  became  a  major  problem  in  PEM  fuel  cells  operation,  the  catalyst 
layer  model  was  modified  by  having  liquid  water  and  polymer  elec¬ 
trolyte  in  its  pores  [17].  The  reactant  gases  were  also  included  in  the 
pores  when  it  became  clear  that  two  phase  flow  may  have  occurred 
in  the  catalyst  layer  [18]  and  in  the  electrode  backing  [19].  The  MEA 
centered  model  was  further  varied  by  incorporating  a  one  dimen¬ 
sional  single  phase  flow  model  of  reactant  gas  in  the  flow  channels 
and  heat  transfer  between  the  bipolar  plate,  cooling  water,  reactant 
gas  and  MEA  [6,7,15].  This  was  followed  by  the  development  of  two 
[20,21]  and  three  dimensional  MEA  models  [22]. 

With  the  advent  of  computational  fluid  dynamic  modelling 
(CFD)  and  its  success  in  modelling  complex  fluid  flow  problems  in 
engineering  at  the  turn  of  the  new  millennium,  water  transport  in 
the  gas  flow  channels  was  modelled  extensively  using  CFD  with  the 
water  and  the  energy  fluxes  at  the  MEA  surfaces  treated  as  simple 
source  terms  at  the  boundary  [23-27].  The  CFD  model  was  further 
extended  by  replacing  the  boundary  water  source  with  Henry’s  Law 
where  the  water  content  in  the  gas  phase  is  a  function  of  the  mois¬ 
ture  content  in  the  electrolyte  [28]  and  by  considering  the  existence 
of  two  phase  flow  in  the  gas  channels  [29-31  ]. 

At  almost  the  same  time,  a  more  unified  approach  to  PEM  fuel 
cells  model  emerged  where  the  electrolyte,  the  catalyst  layer  and 
backing,  and  the  gas  flow  channel  were  rigorously  modelled.  The 
water  and  ion  transport  models  in  the  MEA  were  either  based 
on  the  earlier  MEA  dominated  models  [32]  or  the  more  rigorous 
Stefan-Maxwell  multi-component  diffusion  model  [33,34].  Water 
transport  in  the  gas  channel  was  in  the  gas  phase  only  while  that  in 
the  porous  catalyst  layer  was  thought  to  occur  in  two  homogenous 


phases  [18,32,34]  or  in  multiphase  mixtures  [35-38].  The  water 
transport  in  the  gas  channel  and  the  MEA  were  coupled  by  using 
either  Henry’s  Law  [37]  and  Newton’s  Law  of  cooling  [39]  to  gener¬ 
ate  the  interfacial  water  source  or  using  interfacial  water  flux  from 
the  water  transport  equations  in  the  electrolyte  and  catalyst  layers, 
and  in  the  gas  flowing  in  the  channels  [40].  The  unified  model  was 
further  improved  by  incorporating  two  phase  water  transport  in 
the  channels  [41-44]. 

A  less  known  but  equally  important  PEM  fuel  cells  model,  the 
lumped  or  stirred  tank  reactor  model,  where  the  fuel  cell  is  regarded 
as  a  black  box  and  the  fuel  cell  reaction  is  assumed  uniform  through¬ 
out  the  fuel  cell  without  any  spatial  variation,  was  principally  used 
in  modelling  the  dynamic  response  and  control  of  PEM  fuel  cells 
[45-49].  Its  use  for  water  management  is  limited  because  it  could 
not  show  the  spatial  variation  of  water  content  in  the  PEM  fuel 
cells.  However,  if  the  single  stirred  tank  reactor  or  lumped  model 
were  extended  to  multiple  lumps  or  stirred  tank  reactors  where 
the  fuel  cell  is  divided  into  multiple  stirred  tanks  in  series  along  the 
flow  field,  then  water  content  variation  along  the  channel  could  be 
determined  easily  and  the  accuracy  of  prediction  could  be  adjusted 
by  the  number  of  tanks  used  [50-52]. 

In  this  paper,  a  model  for  the  water  content  profiles  along  the 
flow  channels  at  both  the  anode  and  the  cathode  was  developed 
in  order  to  understand  the  water  flow  in  the  PEM  fuel  cell.  The 
new  model  assumes  that  the  PEM  fuel  cell  is  at  steady  state  and  at 
constant  temperature  and  pressure.  Both  the  air  and  hydrogen  are 
assumed  to  distribute  uniformly  and  to  flow  in  the  laminar  regime. 
The  hydrogen  and  air  velocities  in  the  experiment  were  4  m  s-1  and 
8 ms-1  and  their  Reynolds  numbers  were  785  and  1570,  respec¬ 
tively.  The  flow  regimes  in  both  cases  were  close  to  the  transition 
region  but  still  within  the  laminar  regime.  The  flow  is  also  assumed 
to  be  fully  developed  everywhere  and  the  velocity  is  assumed  con¬ 
stant.  The  pressure  drops  in  both  cases  were  quite  small  at  about 
440  Pa  and  880  Pa,  respectively  because  the  friction  factors  in  both 
cases  were  low  at  0.01  and  0.012,  respectively.  If  hydrogen  depletion 
is  ignored,  the  velocity  would  only  drop  slightly  and  the  constant 
velocity  assumption  holds.  The  hydrogen  stoichiometric  ratio  was 
about  10  and  hydrogen  usage  is  only  about  10%.  Hence  if  hydrogen 
depletion  was  considered,  the  velocity  would  only  drop  by  10%  of 
its  inlet  value.  The  velocity  can  still  be  assumed  constant.  The  gas 
mixtures  are  assumed  to  behave  as  an  ideal  gas.  In  order  to  avoid 
confusion  with  liquid  water,  the  term  moisture  is  used  to  denote 
water  vapour  in  the  gases  and  in  the  membrane.  Since  the  water 
management  strategy  for  flood  prevention  is  by  keeping  the  mois¬ 
ture  contents  of  both  gases  below  saturation,  the  water  transport 
in  the  gas  flow  channel  is  therefore  assumed  to  occur  in  the  gas 
phase  only.  The  single  phase  gas  flow  in  the  channel  can  be  main¬ 
tained  if  the  gas  velocity  is  kept  quite  high  even  when  the  moisture 
in  the  gas  reaches  the  dew  point  [53].  The  moisture  balance  in  an 
infinitesimal  section  of  the  flow  channel  of  width  w,  depth  b  and 
length  dz  at  the  anode  as  shown  in  Fig.  1 ,  assuming  there  is  no  axial 
diffusion,  is  given  by: 


dyA  =  RT , 
dz  vbpjA 


(1) 


where  y A  is  the  molar  moisture  content  in  the  bulk  hydrogen,  JA 
is  the  condensation  flux  of  water  on  the  surface  of  the  anode,  R  is 
the  gas  constant,  T  is  the  temperature,  v  is  the  gas  velocity,  b  is  the 
depth  of  the  flow  channel  and  p  is  the  total  gas  pressure.  The  driving 
force  of  the  moisture  mass  transfer  flux  from  the  hydrogen  gas  to 
the  membrane  by  evaporation  is  given  by  the  difference  between 
yA  and  the  equilibrium  moisture  content  of  the  hydrogen  feed  at 
the  surface  of  the  anode  and  the  anode  temperature,  y* .  As  the 
moisture  content  of  the  membrane  is  low,  it  is  given  by  Henry’s 
law,  y*  =  Hx* ,  where  x £  is  the  equilibrium  moisture  content  of  the 
membrane  [54].  If  water  is  assumed  to  condense  on  the  wet  surface 
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of  the  anode,  then  the  condensation  flux  is  given  by  the  following 
relationship  for  mass  transfer: 


jA  =  kA(yA-y*A)  (2) 

where  kA  is  the  condensation  coefficient.  Eq.  (1)  then  becomes 


dyA 

dz 


RTkA 

vbp 


(yA-yiU 


(3) 


to  evaporate  from  the  wet  surface  of  the  cathode,  then  the  evap¬ 
oration  flux  from  the  wet  surface  of  the  cathode  is  given  by  the 
following  [56,57]: 

Jc  =  kc(ya-yc)n  (7) 

where  yCi  is  the  water  content  of  the  air  near  the  surface  of  the  cath¬ 
ode,  which  is  assumed  to  approximate  Raoult’s  Law  at  the  cathode’s 
surface.  By  virtue  of  Eq.  (5),yci  is  given  by: 


The  analytical  solution  is  given  by 

yA=yA  +  (yAF  -  y*A)  exP  ( ^ )  (4) 

where  Pe  =  vbP/RTkA  is  the  anode  Peclet  number,  z*=z/L  is  the 
dimensionless  distance,  I  is  the  total  length  of  the  flow  field  andyAF 
is  the  hydrogen  feed  moisture  content.  The  membrane  surface  on 
the  anode  side  is  assumed  to  be  saturated  with  water,  but  the  water 
content  at  the  membrane  surface  on  the  cathode  side  is  xci,  which  is 
not  in  equilibrium  with  air.  If  xci  is  smaller  than  x^,  then  there  is  no 
back  diffusion.  At  steady  state,  the  water  content  in  the  membrane 
is  assumed  to  drop  linearly  from  the  anode  to  the  cathode.  This  is  a 
reasonable  assumption  because  the  fuel  cell  was  run  at  low  current 
density  in  such  a  way  as  to  remove  water  formed  by  reaction  at  the 
cathode  as  soon  as  it  was  formed  by  rapid  evaporation  using  low 
humidity  or  dry  air.  The  very  little  accumulation  of  water  on  the 
cathode  side  was  not  sufficient  to  increase  its  concentration  at  the 
cathode  side  of  the  membrane  to  cause  any  significant  back  diffu¬ 
sion  to  the  anode.  In  addition  there  was  a  net  flow  of  water  from 
the  anode  and  cathode.  The  latest  study  on  electro-osmotic  drag 
for  PEM  fuel  cells  operated  at  high  hydrogen  humidity  and  dry  air 
[55]  which  is  similar  to  the  present  study,  confirmed  that  the  mois¬ 
ture  content  dropped  across  the  membrane  from  the  anode  to  the 
cathode  and  the  moisture  profile  was  linear  at  low  current  density. 
Hence,  the  moisture  at  the  surface  of  the  cathode  is  given  by: 


xCi  =  xA  -  M  = 


y*A 

H-M 


(5) 


where  M  is  a  coefficient.  The  driving  force  of  the  moisture  mass 
transfer  flux  from  the  membrane  at  the  cathode  is  given  by  the 
difference  in  the  moisture  content  at  the  surface  of  the  cathode  and 
the  moisture  content  in  the  bulk  air,yA.  Hence,  the  moisture  balance 
in  an  infinitesimal  section  of  the  flow  channel  of  width  w,  depth  b 
and  length  dz  at  the  cathode  as  shown  in  Fig.  1,  assuming  there  is 
no  axial  diffusion,  is  given  by: 


dye 

dz 


RT 

vbP 


C/c  +  r) 


(6) 


where  Jc  is  the  condensation  flux  of  water  on  the  surface  of  the 
anode  and  r  is  the  rate  of  water  produced  by  the  electrochemical 
reaction  at  the  cathode  per  unit  area  of  cathode.  If  water  is  assumed 


dz 
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yci=Kxci=K(^)  (8) 

where  K  is  a  coefficient  and  kc  is  the  evaporation  mass  transfer 
coefficient,  given  by: 

kc  =  D  +  Evm  (9) 


The  values  of  n  and  m  are  either  1.0  and  0.5756,  respectively,  for 
the  Bansal  and  Xie  evaporation  model  [56]  or  0.82  and  1.0,  respec¬ 
tively,  for  the  Tang  and  Etzion  evaporation  model  [57].  The  rate  of 
water  produced  by  the  electrochemical  reaction,  r,  at  the  cathode 
is  given  by: 


(10) 


After  simplifying  the  current  density  expression,  ic,  [58]  it  can 
be  written  as: 

ic  =  ^G(l-yc)  (ID 

where: 

C  =  *4(1  -  v%exp(-K6(</>s-</>m -</&,))) 

X  coth  y/K5  exp(-K6{4>s  -  <£m  -  0£q))  (12) 


„  24te(l  —  g)FDq2 
4= - d2 - 

Up 

_  ioctfsdp 
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Ko=  — 

6  2  RT 


(13) 

(14) 

(15) 


where  te  is  thickness  of  the  GDI,  dp  is  the  pore  diameter,  s  is  the 
dry  porosity,  as  is  the  specific  active  surface  area,  is  the  refer¬ 
ence  molar  concentration  of  oxygen,  </>s  and  0m  are  electrode  and 
membrane  potentials,  respectively,  A0eq  is  the  potential  difference 
between  the  electrode  and  the  membrane,  F  is  the  Faraday’s  con¬ 
stant,  Dq2  is  the  diffusion  coefficient  of  the  dissolved  oxygen  in  the 
electrolyte  and  i0c  is  the  cathodic  exchange  current  density.  Eq.  (6) 
can  be  rewritten  as: 


dye 

dt 


y\  \ 

H-M ) 


+ 


G 

2 Fvb 


(i  -yc) 


(16) 


dye 
dz * 


T(Q-yc)n  +  ^(l-yc) 


(17) 


where  Pe  =  vbP/RTke  is  the  cathode  Peclet  number,  Da  =  G/2Fke  is 
the  Damkohler  number  and  Q.  =  I<  (y*/H  -  m).  If  m  is  small  and 
JO^H,  then  Q.  «  y* .  If  n  and  m  are  not  integers,  then  Eq.  (17)  has  no 
closed  form  analytical  solution  and  can  only  be  solved  by  using  a 
numerical  method  such  as  the  fourth  order  Runge  Kutta  algorithm. 
However,  if  the  adopted  values  n  and  m  are  1.0  and  0.5756  [56], 
respectively,  and  if  the  velocity  is  relatively  constant  throughout 
the  channel,  then  Eq.  (17)  becomes: 


Fig.  1.  Schematic  diagram  of  control  volumes  for  differential  mass  balance  of  water 
in  the  anode  and  cathode  gas  channels  in  a  PEM  fuel  cells. 


dyc  Da  +  Q.  /  Da  +  1  \ 
dz*  =  Pe  V  Pe  )yc 


(18) 
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If  Eq.  (20)  is  integrated  with  respect  toyc»  then: 


yc 


(Da  +  Q) 
(Da  +  1) 


ycF  - 


(Da  +  Q)] 
(Da  +  1). 


exp 


/  Da  + 1 
V  Pe 


) 


z* 


(19) 


Since  the  relative  humidity,  0,  is  the  parameter  that  is  measured 
in  the  experiment,  a  relationship  between  the  mole  fraction  and 
the  relative  humidity  is  needed.  The  mole  fraction  of  water  vapour 
is  given  by: 


y  =  r-f  (20) 

where  ps at  is  the  saturated  vapour  pressure  of  water.  The  moisture 
flux  across  the  membrane,  Jw,  can  also  be  written  in  terms  of  the 
overall  gas  side  mass  transfer  coefficient,  kg: 


Jw  —  kg  (ya  —  yc) 


(21) 


The  overall  gas  side  moisture  mass  transfer  coefficient  kg  for 
both  the  anode  and  cathode  can  be  estimated  from  the  anode  and 
the  cathode  moisture  mass  transfer  coefficients  kA  and  kc: 


1  _  K  1 
kg  ~  kA  +  kc 


(22) 


Both  Eqs.  (4)  and  (19)  were  to  be  validated  by  experimental  data 
from  a  single  cell  PEMFC  that  was  run  at  different  temperatures. 
An  alternative  water  management  strategy  that  forces  the  water  to 
flow  from  the  anode  to  the  cathode  was  implemented  in  the  exper¬ 
iment.  This  was  achieved  by  saturating  the  hydrogen  feed  (0  =  1) 
and  using  dry  air  (0  =  0.005).  Two  pressure  levels  were  used,  1  bar 
and  1.5  bar  at  the  anode  and  the  cathode,  respectively. 


Fig.  2.  Serpentine  flow  field  with  three  holes  for  the  humidity  sensors. 


4.  Materials  and  methods 

4.1  Single  cell  PEMFC 

A  single  cell  PEMFC  was  fabricated  and  used  for  the  model  vali¬ 
dation.  The  specifications  of  the  single  cell  PEMFC  are  given  Table  1. 
Fig.  2  shows  the  graphite  plate  with  a  conventional  serpentine  flow 
field  and  the  positions  of  the  three  humidity  sensors  installed  along 
the  field.  The  positions  of  the  humidity  sensors  on  the  anode  and 
cathode  are  shown  in  Fig.  3.  An  expanded  view  showing  the  various 
components  of  the  single  cell  is  shown  in  Fig.  4.  Two  humidity  sen¬ 
sors  were  also  installed  at  the  feed  points  of  the  air  and  hydrogen 
on  the  single  cell  to  measure  the  external  humidity  of  reactant  gas 
before  entering  the  cell. 

The  temperature  of  the  single  cell  fuel  cell  was  controlled 
by  temperature  controllers  in  the  Fuel  Cell  Test  Station  for  both 
hydrogen  and  air  flows.  The  humidity  sensors  were  read  when  a 
steady  state  of  operation  was  established  after  2  h  of  operation.  A 
schematic  flow  diagram  of  the  experimental  setup  is  shown  in  Fig.  5. 


Table  1 

Single  cell  PEM  fuel  cell  specifications. 


Size 

200  mm  x  200  mm  x  5  mm 

Active  area 

200  cm2 

Active  layer  thickness 

5  fxm 

Specific  surface  area 

1.0  x  107  m2  itt3 

Flow  field  pattern 

Serpentine 

Flow  channel 

Gas  channel  size:  width  =  5  mm, 
depth  =  2.5  mm,  Rib  =  2.5  mm, 
length  =  2989  mm 

Current  collector 

Copper 

End  plate 

Aluminium  =  0.01  m 

Gasket 

Silicone  rubber  0.001  m 

Electrolyte  membrane 

Nafion®  112 

Position  of  humidity  sensor  1 

457  mm 

Position  of  humidity  sensor  2 

1266  mm 

Position  of  humidity  sensor  3 

2532  mm 

4.2.  Experimental  work 

An  Arbin  Fuel  Cell  Test  Station  was  used  to  test  the  single  cell 
PEMFC.  The  temperature  was  varied  from  25  °C  to  60  °C,  while 
the  pressure  at  the  anode  and  the  cathode  was  kept  at  1  bar  and 
1.5  bar,  respectively.  The  flow  rates  of  hydrogen  and  air  were  fixed 
at  3Lmin-1  and  6Lmin-1,  respectively.  The  high  flow  rates  were 
chosen  so  that  single  phase  gas  flow  in  the  channel  was  ensured 
even  when  the  moisture  in  the  gas  reaches  saturation  [53].  Table  2 
shows  the  operational  parameters  of  the  cell.  The  single  cell  fuel 
cell  was  operated  by  externally  humidifying  the  hydrogen  feed  gas 
at  the  anode  until  it  was  saturated  with  moisture  (0  =  1.0).  The  air 
at  the  cathode  was  dried  to  a  low  relative  humidity,  0  =  0.005.  The 
water  profile  models  at  the  anode  and  the  cathode  were  then  vali¬ 
dated  by  the  humidity  profiles  along  the  serpentine  flow  field.  The 
polarization  curves  for  the  single  cell  PEMFC  during  the  experiment 
at  various  temperatures  are  shown  in  Fig.  6. 


Current  collector  Graphite  plates  Insulating  gasket 


Fig.  4.  Expanded  view  of  the  single  cell  PEMFC. 
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Table  2 

Operating  parameters  of  single  cell  PEM  fuel  cell. 


Operating  parameters 

Value 

Cell  pressure 

1.0-1. 5  bar 

Cell  temperature 

25-60 °C 

Hydrogen  flow  rate 

3  Lmin-1 

Air  flow  rate 

6Lmin_1 

Relative  humidity  (anode) 

1.0  RH 

Relative  humidity  (cathode) 

0-0.005  RH 

Exchange  current  density 

10  mA  cm-2 

The  potential  of  electrode 

0.6-0.8  V 

5.  Results  and  discussion 

Fig.  7  shows  the  four-point  moisture  content  profiles  at  four  dif¬ 
ferent  temperatures  along  the  serpentine  flow  field  on  the  anode 
side  that  were  measured  when  a  steady  state  was  reached  after  2  h 


of  operation.  Since  the  hydrogen  feed  was  saturated  with  moisture, 
its  moisture  content  increased  with  temperature.  The  moisture  con¬ 
tent  profile  of  the  hydrogen  decreased  more  rapidly  as  temperature 
was  increased  because  the  mass  flux  rate  of  moisture  through  the 
MEA  increased  more  dramatically  with  increasing  temperature. 

The  similarly  obtained  moisture  content  profiles  on  the  cathode 
side  are  shown  in  Fig.  8.  The  moisture  content  profiles  increased 
more  rapidly  at  a  higher  temperature  because  the  mass  flux  rate  of 
moisture  through  the  MEA  and  the  rate  of  water  production  by  the 
electrochemical  reaction  at  the  cathode  increased  more  drastically 
with  increasing  temperature. 

The  model  equations  were  validated  by  determining  their 
parameters  on  the  anode  and  cathode  sides  using  the  least  squares 
algorithm.  The  experimentally  determined  parameters  are  given  in 
Tables  3  and  4.  The  sums  of  squares  of  residuals  are  all  very  small, 
indicating  that  all  the  experimental  data  closely  fit  the  model  equa¬ 
tions.  These  results  show  that  the  use  of  the  water  evaporation 
model  of  Bansal  and  Xie  [56]  in  the  moisture  profile  model  leads 
to  a  simple  closed  form  solution  that  was  easily  fitted  to  the  exper¬ 
imental  data.  The  evaporation  model  of  Tang  and  Etzion  [57]  did 
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Fig.  6.  Polarization  curves  of  the  MEA  at  constant  temperatures. 


Fig.  7.  Moisture  content  profiles  on  the  anode  side  at  various  temperatures. 
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Table  3 

Parameter  estimation  of  the  model  equation  on  the  anode  side. 


Table  4 

Parameter  estimation  of  model  equation  on  the  cathode  side. 


not  lead  to  a  closed  form  solution  and  was  therefore  not  pursued 
further  in  this  paper. 

The  moisture  profiles  on  the  anode  side  were  influenced  by  the 
Peclet  number,  which  compares  the  relative  strengths  of  convection 
and  mass  transfer  processes.  The  large  value  of  the  Peclet  number 
at  25  °C  compared  to  that  at  higher  temperatures  indicated  that  the 
mass  flux  is  very  small  compared  to  the  convection  of  the  gas  at  this 
temperature.  The  Peclet  number  on  the  anode  side  decreased  from 
59  to  0.32  when  temperature  was  increased  from  25  °C  to  60  °C.  This 
was  caused  not  only  by  the  increasing  moisture  mass  transfer  coef¬ 
ficient  and  moisture  diffusion  in  the  membrane  with  temperature, 
but  also  by  water  molecules  being  dragged  by  the  proton  stream 
across  the  membrane  by  electro-osmotic  drag. 

On  the  other  hand,  the  moisture  profiles  on  the  cathode  side 
were  influenced  by  both  the  Peclet  and  the  Damkohler  numbers. 
The  Damkohler  number  compares  the  relative  strengths  of  the  elec¬ 
trochemical  reaction  and  mass  transfer  processes  at  the  cathode. 
The  Peclet  number  on  the  cathode  side  shows  a  similar  trend  to  the 
Peclet  number  on  the  anode  side.  However,  the  Damkohler  num¬ 
ber  initially  increased  to  a  maximum  at  40  °C,  but  then  fell  back 
to  a  lower  level  at  60  °C.  It  can  be  deduced  that  the  removal  rate 
of  moisture  by  mass  transfer  was  initially  outstripped  by  moisture 
production  from  the  electrochemical  reaction  up  to  40  °C.  Above 
40  °C  the  moisture  removal  rate  by  mass  transfer  overtook  the  rate 
of  moisture  production  from  the  electrochemical  reaction  at  higher 
temperatures  up  to  60  °C. 

Figs.  9-12  show  the  composite  moisture  profiles  at  both  the 
anode  and  the  cathode  at  various  temperatures.  The  moisture  gra¬ 
dient  across  the  MEA  was  initially  large  because  the  hydrogen 
feed  was  saturated  with  moisture  and  the  air  was  dry.  Both  the 
anode  and  the  cathode  moisture  profiles  converged  to  the  same 


asymptotic  value  for  each  temperature.  The  asymptotic  moisture 
content  values  correspond  to  the  equilibrium  moisture  contents  of 
the  membrane  at  the  same  temperature. 

The  anode  moisture  mass  transfer  coefficient,  kA ,  cathode  mois¬ 
ture  mass  transfer  coefficient,  kc,  and  overall  moisture  mass 
transfer  coefficient,  kg,  were  estimated  from  the  Peclet  number  data 
and  are  shown  in  Table  5.  Moisture  fluxes  across  the  MEA  at  various 
temperatures  are  shown  in  Fig.  13. 

The  overall  moisture  mass  transfer  coefficients  and  the  effective 
moisture  diffusivities  are  of  the  same  order  of  magnitude  as  those 
given  in  a  previous  work  [59],  except  at  25  °C.  The  large  discrepancy 
at  25  °C  is  caused  by  the  loss  of  accuracy  of  the  humidity  measure¬ 
ment  at  low  temperature.  The  order  of  magnitude  correspondence 
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Fig.  8.  Moisture  content  profiles  on  the  cathode  side  at  various  temperatures. 
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Table  5 

Estimated  values  of  moisture  mass  transfer  coefficient  and  effective  diffusivity. 


kA  (mol  cm-2  s_1) 

kc  (mol  cm-2  s_1) 

K-values  for 
moisture  [16] 

Experimental  kg 
(mol  cm-2  s-1) 

Published  kg 
(mol  cm-2  s_1)  [22] 

Experimental  DE 
(mol  cm-2  s_1) 

4.66  x  10“5 

2.46  x  10-5 

54.93 

8.19  x  10-7 

1.35  x  10“5 

1.64  xlO-8 

2.88  x  10“3 

2.97  x  10“3 

64.43 

4.41  x  10“5 

2.85  x  10“5 

8.81  x  10“7 

3.56  x  10“3 

3.00  x  10“3 

70.75 

4.95  x  10“5 

3.60  x  10“5 

9.99  x  10“7 

7.60  x  10“3 

2.94  x  10“3 

77.08 

9.54  x  10“5 

4.70  x  10“5 

1.91  x  10“6 

Fig.  12.  Moisture  content  profiles  in  the  anode  and  cathode  at  60  °C. 

of  the  overall  mass  transfer  coefficients  and  the  effective  moisture 
diffusivities  with  published  data  validates  the  model. 

Moisture  flux  increased  with  increasing  temperature,  but 
decreased  along  the  flow  field  and  vanished  at  the  outlet  of  the  flow 


field  when  the  moisture  content  reached  the  equilibrium  moisture 
content  of  the  membrane  at  the  outlet.  Water  formed  by  the  reac¬ 
tion  and  moisture  from  the  anode  was  immediately  removed  by 
the  dry  air  because  of  the  large  moisture  content  gradient  between 
the  cathode  and  the  dry  air.  The  moisture  did  not  reach  its  dew 
point  at  any  point  in  the  cathode,  and  no  water  droplets  that  may 
have  caused  flooding  were  formed  in  the  cathode.  The  conventional 
strategy  of  saturating  both  the  hydrogen  and  air  feeds  was  found  to 
cause  rapid  formation  of  water  droplets  in  the  cathode  side,  flood¬ 
ing  and  blocking  the  flow  field  channels  [60]  and  the  electrode  [61  ], 
because  moisture  dew  point  was  already  reached  near  the  air  inlet 
[61]. 

6.  Conclusions 

The  moisture  profile  in  the  flow  field  channels  of  the  anode 
was  shown  to  depend  on  the  moisture  Peclet  number,  the  decrease 
of  which  with  increasing  temperature  indicates  increasing  mois¬ 
ture  mass  transfer.  On  the  other  hand,  the  moisture  profile  in  the 
cathode  was  shown  to  depend  on  both  the  Peclet  number  and  the 
Damkohler  number.  The  trend  of  the  Peclet  number  in  the  cath¬ 
ode  followed  closely  that  of  the  anode.  The  Damkohler  number 
decreased  with  temperature,  indicating  increasing  moisture  mass 
transfer.  The  moisture  profile  models  were  successfully  validated 
by  the  published  data  of  the  estimated  overall  mass  transfer  coef¬ 
ficient  and  moisture  effective  diffusivity  with  those  in  published 
data  of  the  same  order  of  magnitude.  The  strategy  of  saturating  the 
hydrogen  feed  and  using  dry  air  in  the  present  work  has  been  shown 
to  successfully  prevent  water  droplet  formation  in  the  cathode  and 
hence  prevent  flooding. 
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